N=32;
j=sqrt(-1);
a=ones(N,1);
row=0;
for i1=1:2:N;
	a(i1) = sqrt(row+i1-1+.77777);	
	a(i1+1) = (row+i1-1+.22222)*(row+i1-1+.22222) / N - N/2;	
end;
a(1)=N+3;
a(3)=1-N;
A=fft(a);
A=A(1:N/2);
if N<=32
	A
end
max(abs(A-atrans)./(max(abs(A))/N))
